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Abstract 


Host  studies  of  Si  oxidation  comasnca  with  a  discussion  of 
ths  Linear -Parabolic  oxidation  model  developed  by  a  number  of 
workers  in  the  1960 's.  The  limits  of  the  model  are  pure  diffusion 
of  oxidant  for  thick  SiO*  films  and  a  surface  reaction  limitation 
for  thin  films.  The  steady  state  picture  of  this  series  reaction 
scheme  is  discussed  and  used  to  explain  new  experimental  results. 
New  data  relevant  to  Si  oxidation  is  presented  on  the  following 
subjects!  five  orientations  of  Si;  photonic  excitement;  intrinsic 
film  stress;  sillcide  oxidation.  The  role  of  electrons  on  the 
oxidation  kinetics  is  elucidated.  A  thermionic  emission  model  for 
the  initial  stages  of  oxidation  is  proposed. 


Introduction 


% 
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The  linear  parabolic,  L-P,  oxidation  model (1-4)  is  typically 
adopted  as  the  basis  from  which  to  describe  Si  oxidation 
phenomena.  This  model  is  based  on  a  consideration  of  mass  balance 
and  the  arrival  at  the  steady  state  of  essentially  two  rate 
processes.  One  of  the  processes  is  the  transport  of  oxidant 
across  the  SiO*  film,  represented  in  Figure  1  as  Fa ,  and  the 
other  process  is  the  reaction  of  oxidant  with  Si  at  the  Si-SiO* 
interface  which  is  shown  as  F* .  It  is  intuitive  that  these 
processes  are  essential  for  oxidation.  Which  of  these  processes 
is  kinetically  dominant  is  not  intuitive  and  the  subject  of 
numerous  experimental  investigations.  Before  specific 
experimental  studies  are  considered  which  elucidate  this  issue, 
it  is  useful  to  consider  the  implications  of  mass  balance  and  the 
steady  state. 

Diffusion  of  Oxidant  in  SIO. 

Fickian  diffusion  is  usually  taken  to  mean  that  the  mass 
flux,  F« ,  is  proportional  to  a  concentration  gradient,  dC/dx, 
which  in  the  steady  state  is  time  independent: 

F,  »  DdC/dx  (1) 

where  D,  the  diffusivlty,  is  the  constant  of  proportionality.  A 
number  of  concordant  studies(5,6)  show  that  the  0*  flux  through 


1 


SiOa  is  overwhelmingly  without  interaction  with  the  SiO*  network. 
On  the  other  hand,  for  the  case  where  even  traces  of  H«0  are 
present  the  exchange  of  0  with  the  network  is  rapid <6, 7)  and  the 
oxidation  kinetics  are  decidedly  altered ( 8, 9 ) .  The  overall 
oxidation  rate  has  been  shown  to  scale  with  oxidant 
pressure ( 1®,  11 ) .  These  kinds  of  oxidation  studies  along  with  the 
intuitive  idea  that  as  the  oxide  grows  the  reactants  are 
increasingly  separated,  yield  the  conclusion  that  the  transport 
of  oxidant,  or  the  oxidant  supply  function  is  kinetically 
important. 

The  Interf aqlal  flffgUgn 

It  is  clear  that  SiOt  cannot  form  without  the  reaction  of  Si 
and  0  atoms.  The  reaction  can  be  expressed  in  terms  of  a  flux, 

F« ,  where  the  concentrations  of  both  chemical  species  are 
important : 

F,  >  kCOHSil  (2) 

and  where  k  is  a  rate  constant.  Since  for  a  specified 
crystallographic  orientation,  the  availability  of  Si  atoms  on  the 
surface  is  constant,  [Si]  could  be  included  in  k.  However,  it  is 
now  well  established  that  the  oxidation  rate  is  dependent  upon 
the  orientation  up  to  thousands  of  angstroms < 12 )  and  thus  it  is 
explicit  in  equation  <2).  We  ignore  that  CO]  and  CSi]  may  appear 


in  equation  <2>  raised  to  powers  other  than  unity  and  that  oxygen 


molecules,  Oa »  may  ba  important. 

Series  Procesaea-Masa  Balance-Steady  State 

With  the  processes  of  transport  of  oxidant  and  subsequent 
reaction  is  series,  mass  balance  requires  that  the  fluxes  relax 
to  the  same  value,  F,  as: 


F  =  F,  =  F,  (3) 

This  relaxation  will  cause  the  establishment  of  a  quasi 
stationary  concentration  at  the  Si-SiO*  interface,  C« ,  as: 

dC./dt  *  0  (4) 

This  physical  picture  is  shown  in  Figure  2  after  Levich(44)  where 
F(C* )  for  both  F,  and  F«  is  graphed  versus  C« .  It  is  seen  that  at 
the  0« -SiO,  interface,  the  maximum  value  for  F,  is  calculated  at 
Ca  =  0  as : 


F,  (0)  =  DC,  /L  (5) 


where  C,  is  the  solubility  of  oxidant  in  the  SiO«  and  L  is  the 
SiO*  thickness.  F,  goes  to  0  at  C«  =  C,  and  varies  linearly  in 
between,  hence  a  straight  line  is  drawn  between  these  extrema  in 
F, .  F«  can  have  any  maximum  value  and  is  shown  with  both  a  larger 
and  smaller  maximum  value  than  F, .  The  maximum  value  for  F«  is  at 
C»  *  C,  and  has  a  value  of  0  at  C«  =  0.  These  extrema  in  F(  can 
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be  connected  by  a  line.  The  intersection  of  the  two  functions 


yields  the  steady  state  solution  for  C* ,  viz.  the  C«  value  that 
satisfies  equation  <3>  above. 

While  in  general  the  maximum  value  for  F»  can  be  smaller  or 
larger  than  Ft ,  below  we  show  evidence  which  supports  the  former. 
Also,  it  is  seen  that  the  F  value  obtained  from  the  steady  state 
must  be  less  than  the  smallest  flux,  since  the  intersection  for 
two  downwardly  sloping  lines  yields  such  an  intersection.  The 
slope  of  the  F*  function  is  determined  by  k  for  any  value  of  Ct , 
and  the  slope  of  Ft  by  D.  Changes  in  oxidation  rate  can  be 
analyzed  by  noting  the  changes  in  the  intersection  point  as  a 
result  of  changes  in  the  parameters  in  each  flux  expression.  We 
will  use  this  idea  later  to  explain  some  experimental  results. 
Thus  Figures  1  and  2  provide  pictures  of  the  important  overall 
processes  and  how  these  processes  interact  in  series. 

Results 


The  Dominant  Rate  Process 

In  order  to  address  the  question  posed  above  as  to  which 
rate  process  is  dominant,  several  simple  calculations  are 
relevant. 

First,  a  comparison  is  made  of  the  maximum  0«  diffusive  flux 
(at  C*  ■  0),  Ft  (0),  with  the  experimentally  obtained  oxidation 
rate  converted  to  an  flux,  F(exp).  This  comparison  is  easily 
made  by  forming  a  ratio,  R,  as: 
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R  =  F, <0> /F(exp)  (£) 

Table  1  shows  some  recently  reported  values ( 12, 13 ) .  Literature 
values  for  D(14>  and  Ct (3)  were  used.  It  is  seen  that  for  any 
reasonable  oxide  thicknesses,  R  values  considerably  greater  than 
unity  are  found.  R  values  of  unity  indicate  that  the  reaction 
scheme  is  under  transport  control  which  is  approached  for  large 
SiOt  thicknesses  greater  than  2000  A.  For  R  greater  than  unity, 
ie.  most  of  the  important  SiO«  film  thicknesses,  F  <  Ft  (0)  and 
this  indicates  that  a  steady  state  value  for  Ct  is  non  zero. 

Using  the  D  values  for  □• ,  a  random  walk  calculation  shows 
that  the  net  distance  travelled  by  an  0«  molecule  in  SiOa  in  one 
second  is  quite  large  and  as  shown  in  Table  2  is  about  10*  A. 

Thus  Ft  appears  to  be  quite  large.  A  recent  report  (15)  has  shown 
that  the  thermodynamic  potential  associated  with  oxidant 
transport  during  Si  oxidation  is  smaller  than  the  reaction 
potential  by  more  than  a  factor  of  50. 

Now  considering  these  three  calculations  and  the  picture  in 
Figure  2,  one  can  conclude  that  the  observed  oxidation  rate  as  F 
is  nearer  to  the  F*  line.  This  indicates  a  large  slope  for  the  Ft 
line  (a  large  D)  and  that  Ft  >  F«  for  the  maximum  values. 
Experimental  support  for  this  is  now  presented. 

The  Effect  of  Si  Orientation 

It  has  long  been  recognized  that  the  different  Si 
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orientations  oxidize  at  differing  oxidation  rates<2, 3, 16-18).  For 
the  three  lowest  index  major  Si  orientations,  the  reported  rate 
order  is  as  follows: 

Rate<110)  >  Rate(lll)  >  Rate(100) 

for  steam  at  1  atm<2)  and 

Rate (111)  >  Rate (110)  >  Rate(100) 

for  high  pressure  steam  oxidation ( 16 ) .  Such  differences  were 
accounted  for  in  the  L-P  model  by  adjustments  in  the  linear  rate 
constant,  k,  in  equation  (2).  All  the  authors  agree  that  the 
number  of  Si  atoms  per  area,  Na  «  is  the  important  parameter  and 
this  is  shown  explicitly  in  equation  (2).  Na  t  was  calculated  for 
the  three  major  orientations < 16 )  and  found  to  be  in  the  order: 

Na i  (110)  >  Na  t  (111)  >  Na  t  ( 100 ) . 

Thus  for  one  of  the  above  cited  studies  (high  pressure  steam) 
another  parameter  was  required  to  be  invoked  that  could  dominate 
the  kinetics.  This  was  chosen  to  be  steric  factors  and  a 
geometric  argument  was  proposed  to  explain  the  experimental 
results ( 16 ) . 

Recent  detailed  kinetics  studiest 19-21 )  using  the  three 
major  Si  orientations,  have  conclusively  shown  that  for  1  atm  dry 


□a  oxidation*  th*  initial  rate  order  (leas  than  several  hundred 


A)  is: 


Rate (  110  )  >  Rate  (111)  >  Rate(100> 


which  changes  for  thicker  films  to  the  order: 


Rate(lll)  >  Rate  (110)  Rate(100> 


A  model  was  proposed  for  the  initial  order  and  the  'crossover 


ef f ect " ( 22 ) .  For  the  initial  regime  the  rate  simply  scales  with 


Na  4  as  given  above,  since  the  very  initial  regime  ought  to  be 


dominated  by  the  interface  reaction  where  : 


[Si]  oc  N, 


The  crossover,  which  occurs  at  greater  film  thicknesses  was 


attributed  to  a  film  stress  effect  that  is  manifest  as  the  force 


which  grows  with  the  film  thickness.  The  intrinsic  stress 


develops  as  a  result  of  the  change  in  molar  volume,  AV,  which 


occurs  from  the  conversion  of  Si  to  SiO*  and  is  about  120%.  The 


orientation  dependence  was  thought  to  arise  from  the  application 


of  Hookes  law  as: 


c  =  Ee  (7) 


VJf*1 


vhtrt  o  la  tha  stress,  and  e  i a  the  strain,  and  is  proportional 


to  AV,  and  E  is  Youngs  modulus  which  is  orientation  dependent. 

The  variation  of  E  is: 

E(lll)  >  E ( 110)  >  E(100> 

thus  paralleling  the  order  of  the  oxidation  rate,  R,  above  the 
crossover.  The  stress  is  compressive  in  the  oxide,  but  tensile  in 
the  Si  surface.  This  latter  stretched  bond  situation  in  the  Si 
was  considered  to  be  responsible  for  the  enhanced  R  for  the  (111) 
orientation  above  the  crossover.  It  should  be  noted  that  at  the 
time  this  stress  argument  was  developed,  the  experimental 
measurements  of  SiQa  stresses  on  various  Si  orientations  were  not 
available.  Recently,  a  large  number  of  precise  SiOa  film  stress 
measurements  have  been  reported ( 23, 24 >  which  definitely  show 
that  the  intrinsic  stress  does  not  follow  the  order  for  E  above. 
In  fact.  Figure  3  shows  that  the  measured  stress  on  the  (111) 
surface  is  lower  than  for  the  other  three  orientations  which  are 
nearly  equal.  Thus  the  stress  based  model  as  proposed  above 
cannot  be  correct.  The  fact  that  the  (111)  measured  stress  is  the 
lowest  with  the  stresses  for  the  other  orientations  grouped, 
suggests  an  anomaly  with  only  the  (111)  surface.  The  smaller 
compressive  oxide  stress  for  SiOa  on  the  (111)  surface  may  result 
in  a  larger  D  value  for  0a  thereby  explaining  the  increased 
oxidation  rate.  It  is  this  idea  that  is  used  later  to  explain  the 
orientation  results.  Before  presenting  that  argument  it  is  useful 


to  present  further  more  complete  and  very  recent  oxidation  data 
as  a  function  of  the  Si  substrate  orientation ( 12 ) . 

The  oxidation  data  of  Haseoud  et  al( 19-21)  has  recently  been 
extended  down  to  600°C(12)  where  the  intrinsic  stresses  are 
larger ( 22-26 )  and  now  includes  the  (311)  and  (511)  orientations 
as  well  as  the  (111),  (110),  and  (100)  major  Si  orientations.  The 

crossover  in  the  thickness  time  plots  betweeen  the  (110)  and 
(111)  orientations  has  been  observed  for  all  oxidation 
temperatures  between  1100°C  and  750°C  inclusive.  For  lower 
oxidation  temperatures,  the  films  were  not  grown  thick  enough  to 
exhibit  a  crossover  in  the  thickness  time  plots,  but  a  crossover 
in  the  oxidation  rates  (the  slope  of  the  thickness-time  plots) 
was  observed.  The  (311)  and  <511)  orientations  were  originally 
added  to  the  study  because  they  were  thought  to  have  a  higher 
areal  density  of  Si  atoms,  if  one  assumes  that  the  planes  are 
flat.  However,  this  assumption  was  found  to  be  erroneous ( 27 ) .  The 
(311)  and  (511)  are  both  vicinal  planes  of  the  (100)  and  (111) 
planes.  This  is  shown  in  two  dimensions  to  simplify  the  situation 
with  the  aid  of  Figure  4  in  which  the  angle  for  the  (311)  and 
(511)  planes  are  made  using  combinations  of  the  (100)  terraces 
and  (111)  risers.  Using  the  appropriate  ratios  of  (100)  and  (111) 
planes  the  actual  areal  densities  of  Si  atoms  on  the  various 
surfaces  are  calculated ( 12, 2S )  and  shown  in  Table  3.  The 
oxidation  data  at  700°C  and  1000°C  are  shown  in  Figure  5.  A 
crossover  in  only  the  oxidation  rate  is  obtained  for  the  lower 
temperature,  since  the  data  have  not  been  carried  out  to 


sufficient  SiO«  thickn 


to  obatrvc  the  crossover  on  the 


thickness-time  plot.  For  the  1 000°C  data  the  crossover  is  noticed 
on  the  plot.  For  the  three  major  orientations  the  crossover  has 
been  reported < 19-22 )  and  is  shown  in  Figure  6  for  an  800®C 
oxidation.  In  the  pre  crossover  regime  the  rate  order  is: 

Rate (110)  >  Rate (111)  >  Rate<311>  >  Rate(511>  >  Rate(100) 

This  order  parallels  the  order  for  the  number  of  Si  atoms  as 
given  in  Table  3,  but  agreement  is  not  quantitative.  Table  4 
shows  a  comparison  for  the  700®C  and  1000*0  data  and  while 
agreement  is  better  at  the  higher  temperatures  quantitative 
agreement  with  N« ,  is  lacking.  At  the  present  time  a  proven  model 
that  both  explains  the  crossover  and  is  quantitative  is  not 
available.  However,  a  proposal  will  be  offered  which  appears 
consistent  with  the  available  data  and  the  steady  state  ideas 
presented  above. 

For  an  explanation  of  the  crossover,  we  return  to  a  steady 
state  diagram  shown  in  Figure  7  which  is  similar  to  Figure  2  with 
the  addition  of  some  new  flux  lines.  On  the  right  of  the  diagram 
the  two  different  lines  represent  different  interfacial  fluxes 
due  to  the  different  N« ,  for  the  (110)  and  (ill)  Si  surfaces. 

From  equation  <2),  kCSi]  is  different  according  to  the  N» , .  On 
the  left  is  an  additional  flux  line  representing  the  effect  of 
Increased  compressive  stress  for  all  the  orientations  except  the 
(111)  which  shows  the  smallest  stress.  A  compressive  film  stress 
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decreases  0  for  the  film,  thereby  decreasing  the  F  intercept  and 
the  slope.  The  effect  is  a  decrease  in  the  observed  flux,  ie.  the 
oxidation  rate,  as  indicated  on  the  figure  by  a  flux  line 
intersection  at  a  lower  F  for  the  (110)  plane  relative  to  the 
(ill)  Si  surface.  The  stress  effect  on  the  oxidation  rate 
increases  as  the  film  grows,  because  the  resultant  force 
increases  with  film  thickness.  Since  the  effect  of  stress  on  the 
oxidation  flux  or  D  is  not  yet  quantitative,  we  only  obtain  the 
direction  for  the  change.  Thus  the  reduced  stress  on  the  (111)  or 
the  increased  stress  on  the  other  orientations  causes  a 
corresponding  increase  in  0  for  the  (111)  or  decrease  for  the 
other  orientations.  This  stress  change  causes  a  change  in  the 
steady  state  position  and  hence  the  observed  oxidation  rate 
crossover. 

In  summary  then  the  steady  state  analysis  can  yield  a 
reasonable  if  not  yet  quantitative  idea  of  the  various  oxidation 
parameters. 

A  number  of  recent  studies  have  pointed  out  that  both 
thermal  and  photonic  effects  are  at  play  when  oxidation  is 
performed  in  the  presence  of  high  intensity  light ( 29-34 ) .  There 
exists  evidence  that  photon  energies  greater  than  3  eV  yield 
large  oxidation  rate  enhancements ( 30, 34 ) .  From  the  energy  band 
picture  for  the  Si-SiO*  interface,  as  seen  in  Figure  8,  about  3 
eV  is  the  energy  distance  from  the  Si  conduction  band  to  the  SiO« 
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conduction  band.  Further  oxidation  rate  enhancements  were 
observed  for  photon  energies  greater  than  5  eV  which  corresponds 
to  the  dissociation  energy  for  the  molecule.  It  thus  appears 
that  electron  excitation  plays  a  significant  role  in  the 
oxidation  of  Si,  in  addition  to  the  chemical  effects  of  CO]  and 
CSil  in  the  traditional  oxidation  model. 

Sillclde  Oxidation  Studies 

The  role  of  electrons  on  the  Si  surface  is  further 
substantiated  from  recent  studies  of  the  oxidation  behavior  of  a 
variety  of  metal  silicides.  Firstly,  new  oxidation  data  are 
presented.  A  variety  of  metal  silicides  were  cleaned (35)  and 
oxidized  so  as  to  produce  only  SiO«  and  preserve  the  silicide. 
This  can  be  accomplished  in  the  appropriate  temperature  regime 
(about  800°C)  and  most  importantly  with  the  silicides  deposited 
upon  a  Si  substrate  which  insures  a  supply  of  Si  for  the 
oxidation  reaction,  such  that  the  silicide  can  be  preserved.  Thus 
there  is  a  net  flux  of  Si  from  the  substrate  to  the  silicide-SiOc 
interface  and  a  net  flux  of  metal  towards  the  Si.  The  actual  flux 
may  be  more  complicated  and  found  to  be  different  for  different 
silicides( 36-39 ) .  The  results  of  dry  oxidations  at  800°C  are 
shown  in  Figure  9(40).  Three  bands  of  oxidation  behavior  are 
characteristic:  the  fastest  oxidizing  silicides  are  the 
transition  metal  silicides;  the  intermediate  oxidizing  silicides 
are  the  refractory  silicides;  the  slowest  are  the  semiconducting 
silicides  which  display  oxidation  rates  not  much  different  from 
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SI  itself.  It  was  also  found  that  the  free  carrier  concentrations 
and  the  optical  absorbtion  index,  k,  which  is  related  to  the  free 
carrier  concentration,  all  display  the  same  order  as  the 
oxidation  rate.  In  addition  is  a  report (41 >  that  Si  surfaces  that 
have  a  comparatively  large  density  of  electron  states  near  the 
Fermi  level  oxidize  comparitively  rapidly.  These  experiments  were 
done  under  UHV  conditions  using  disordered  (high  density  of 
states)  and  ordered  (low  density  of  electron  states)  noble  metal 
sub-monolayers  on  single  crystal  Si  substrates.  Taken  together 
these  experimental  findings  point  to  a  profound  oxidation  rate 
enhancment  effect  due  to  the  high  electron  densities.  This  is  in 
further  corroboration  to  the  photonic  results  above. 

A  Thermionic  Emission  Model  For  Si  Oxidation 

Referring  back  to  Figure  8,  arid  now  considering  the 
importance  of  electron  availability,  a  straightf oward  dependence 
of  the  availability  of  electrons,  as  the  electron  flux,  F. , ,  on 
the  oxidation  rate  can  be  tested (42)  using  the  thermionic 
emission  relationship,  viz.  the  Richardson-Dushman  equation: 

F.  ,  =  AT*  exp  (  /kT  )  (8) 

where  A  ia  the  Richardson  constant,  T  is  absolute  temperature,  x* 
is  the  barrier  height  which,  as  typically  applied  to  metals, 
extends  from  the  Fermi  level  of  the  metal  to  vacuum.  For  the 
present  case  of  electron  emission  from  Si  to  SiO# ,  we  consider 
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the  highest  levels  in  Si  that  have  sufficient  electrons  to  supply 
the  oxidation  reaction.  Calculations  have  shown  that  there  are 
more  than  sufficient  0 a  molecules  on  the  SiOa  side  of  the 
interface  as  a  result  of  a  fast  diffusive  flux,  and  more  than 
sufficient  electrons  available  on  the  Si  side (42).  Thus  the  rate 
limiting  step  is  assumed  to  be  emission  of  electrons  over  the 
barrier,  step  2  below,  and  the  entire  oxidation  reaction  is 
envisioned  as  follows: 

1.  e*  (C.B.  of  Si)  -»  e*  (Si  surface) 

2.  e*  (Si  surface)  ■*  e*  (free  in  SiOa  ) 

3.  e-  (in  SiOa  )  *  □«  (from  diffusion)  •+  0«  * 

4.  *  •*  0  ♦  O' 

5.  Si*  (Si  surface)  ♦  0  ♦  0*  ♦  -►  SiO, 

Steps  1,  3-5  are  relatively  fast  with  respect  to  step  2.  Step  1 
represents  the  flux  of  electrons  to  the  Si  surface.  Step  3 
represents  the  thermodynamically  favored  process  of  electron 
attachment  to  0a(43).  This  step  destabilizes  0a  and  leads  to 
decomposition  to  a  neutral  and  charged  0  species (43).  Then  step  5 
considers  the  rapid  reaction  of  the  atomic  and  charged  species. 
Notice  that  Si  is  represented  as  Si*,  since  electrons  were  lost 
from  the  surface. 

One  way  to  test  the  rate  limitation  of  Step  2  is  to  equate 
the  thermionic  emission  electron  flux,  F. . ,  with  the  observed 
oxidation  flux,  F(exp),  as  an  SiOa  flux,  and  then  calculate  the 


value  of  x*  ncctaaary  for  the  equivalency,  and  compare  this  with 


the  energy  band  diagram  and  known  values  for  the  various 


barriers.  Table  5  reports  these  x«  values  obtained  from  oxidation 


data  at  various  temperatures  and  for  the  Si  (110)  and  (100) 


orientations.  These  orientations  represent  the  fastest  and 


slowest  oxidizing  surfaces  in  the  initial  regime.  Also  while  the 


initial  oxidation  rates  were  emphasized,  since  thermionic 


emission  is  expected  to  control  the  very  initial  rate,  larger 


thickni 


rre  included  to  show  trends.  It  is  seen  that  within 


a  few  tenths  of  an  eV  the  value  of  3  eV  is  obtained  for  x»  which 


corresponds  to  the  barrier  between  the  Si  conduction  band  and 


SiOa ,  thus  confirming  the  models  consistency.  It  is  easily 


verified  that  the  Si  conduction  band  contains  sufficient 


electrons  to  supply  the  oxidation  reaction  at  any  temperatures 


above  100°C.  At  room  temperature  it  is  universally  observed  that 


a  thin  native  oxide  film  forms  on  freshly  etched  or  cleaved  Si 


surface  with  a  thickness  that  may  range  from  5  to  20  A. 


Thermionic  electrons  do  not  seem  available  for  this  process. 


However,  when  the  a.10‘  *  cm**  surface  electron  states  on  the  Si 


surface  are  considered,  an  explanation  emerges.  If  each  of  these 


states  can  contribute  a  trapped  electron  to  oxidation,  as  the 


state  is  removed  via  oxidation,  then  this  number  of  electrons 


would  explain  the  formation  of  about  10  A  SiO(,  considering  one 


electron  per  SiO«  molecule  as  per  the  reaction  scheme  above  via 


electron  attachment  to  0. .  This  mechanism  would  cease  once  the 


states  are  destroyed,  as  is  observed  for  the  formation  of  the 


native  oxide.  Then  the  thermionic  supply  function  is  invoked  to 
continue  the  oxidation  at  higher  temperatures.  After  the  initial 
regime  the  reaction  mechanism  may  be  controlled  by  other 
processes  which  would  shift  the  steady  state  point  as  discussed 
above.  A  thickness  dependent  thermionic  emission  does  not  seem 
plausible. 


Summary  and  Conclusions 

Several  new  ideas  and  experimental  results  were  discussed 
pertaining  to  Si  oxidation,  especially  in  the  initial  oxidation 
regime.  Many  of  these  ideas  have  only  recently  or  not  yet 
appeared  in  the  literature. 

One  recent  study  clearly  linked  the  Si  oxidation  rate  to  the 
areal  density  of  Si  atoms  on  the  Si  surface,  N« , .  This  dependence 
existed  to  surprisingly  large  oxide  thicknesses  of  several 
thousands  of  A.  Using  recent  intrinsic  SiO*  stress  measurements, 
the  oxidation  rate  crossover  phenomena  between  the  (110)  and 
(111)  orientations  was  explained.  This  was  understood  using 
pictorial  representations  of  the  essential  fluxes  for  the 
oxidation  process,  and  the  steady  state  solution  for  the  fluxes 
in  a  series  process.  New  results  for  the  oxidation  of  metal 
silicides  to  produce  SiO*  films  along  with  recently  reported 
photonic  effects  on  Si  oxidation  have  pointed  to  the  importance 
of  electron  supply  to  the  Si  oxidation  mechanism.  A  simple 


analysis  of  the  initial  oxidation  rates  in  terms  of  a  thermionic 


•mission  model  yielded  correct  values  for  the  Si-SiOa  barrier. 

Quantitative  aspects  of  several  of  the  presented  ideas  are 
still  lacking  and  require  more  research.  The  area  of  potential 
electronic  implications  of  these  new  studies  has  not  been 
adressed  in  the  paper,  but  such  relationships  likely  exist.  For 
example,  the  origin  of  interface  fixed  charge  may  be  explained  by 
the  emission  of  electrons  into  SiQa  during  oxidation.  Continuing 
research  at  our  laboratory  focuses  on  quantification  of  the 
results  and  the  experimental  verification  of  possible 
relationships  between  materials  parameters  and  electronics 
properties  of  the  Si-SiOa  interface. 
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Figure  1.  The  essential  rate  processes  for  Si  oxidation:  F,  is 

the  flux  of  oxidant;  and  F«  the  reaction  between  Si  and 
oxidant. 

Figure  2.  A  flux  diagram  shoving  the  variation  of  Ft  and  F«  with 
C« ,  the  concentration  of  oxidant  at  the  Sl-SlO* 
interface. 

Figure  3.  Intrinsic  SiQ«  film  stress  as  a  function  of  oxidation 

temperature  and  Si  substrate  orientation  (ref.  24,  with 
permission  of  the  Am.  Vac.  Soc.  ). 

Figure  4.  A  two  dimensional  view  of  (311)  plane  relative  to 

(100)  and  (111)  displaying  the  stepped  nature  of  the 
(311)  surface. 

Figure  5.  Oxidation  data  for  five  Si  orientations  at  a)  700°C  and 
b)  1000°C  (ref.  12  with  permission  of  the  Electrochem. 
Soc.  > . 

Figure  6.  Oxidation  data  at  800°C  for  three  major  orientations, 
<(110)  triangles,  (111)  squares,  (100)  circles)  a)  a 
close  up  of  the  initial  data  prior  to  crossover  and  b) 
after  the  crossover  (ref.  22  with  permission  of 
Electrochem.  Soc.  >. 

Figure  7.  Flux  versus  Ct  for  the  (110)  and  (111)  Si 

orientations,  F# 's,  and  for  two  different  0  values, 

F,  'a. 

Figure  8.  Energy  band  diagram  for  Si-SiO*  interface  shoving  the 
barrier  of  about  3. 15eV  from  the  bottom  of  the  Si 
conduction  band. 

Figure  9.  Oxidation  data  for  a  variety  of  metal  silicides  (ref. 
40). 
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